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STRESSES AND DISPLACEMENTS IN THIN SHELLS 
COMPOSED OF CYLINDRICAL AND SPHERICAL SEGMENTS 
by 
Mario G. Salvadori, A.M. ASCE? 
SYNOPSIS 

Membrane stresses and displacements due to uniform load 
are derived for thin shells in the shape of a half-barrel 
closed at the ends by two quarter-spheres. 

Bending stresses are determined at the intersection of 
the cylinder with the sphere and at the horizontal boundary of 
the composite shell. 

The results are easily extenddd to other types of shells 
built by segments of cylinders and spheres, to dead load 


stresses and to thermal stresses. 


1. INTRODUCTION 

Thin shells of reinforced concrete have been used with 
rapidly increasing popularity to build domes and barrel roofs. 
In recent years both in the United States and abroad new shell 
shapes have been utilized in modern architectural design and 
the analysis of thin shells has become an important chapter 
of structural analysis. 

Barrel roofs are usually supported by transverse stiff- 
eners in the form of arches. The present analysis deals, in- 


stead, with barrel roofs stiffened by quarter-spheres at the 


ends, and is mainly concerned with the interaction between the 


cylindrical and the spherical shell segments under uniform load. 
The membrane stresses and corresponding. displacements sre first 
derived and the bending stresses are then evaluated at the 
intersection of cylinder and sphere. These stresses are of a 
local character and affect only a small portion of the shell. 

] Professor of Civil Engineering, Columbia University, 
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Methods for the evaluation of the bending stresses due to 
the boundary conditions are available and are discussed in de- 
tail. 

The analysis may be easily extended to other types of 


loading 


conditions, among them dead load and temperature 
changes, and to different types of shells made of spherical 


end cylindrical segments. 


2. CYLINDRICAL SHELL UNDER UNIFORM LOAD 


The circular cylindrical shell of Fig. 1 carries a uniform 


vertical load. We let: 
Ny transverse force per unit length, positive when tensile 


= longitudinal force per unit length, positive when tensile 


shear force per unit length, positive in the direction 
of positive x and ¢ 


uniform load per unit of norizontal area 
radius of cylinder 

length of cylinder 

angle measured from crown 


longitudinal coordinate measured from middle of 
cylinder 


= x/L = non-dimensional] longitudinal coordinate 
Under the assumption that the shell thickness be so small 
as to make the flexural rigidity of the shell negligible, 
load p is carried by membrane stresses given by: 
Ny =-pa (L/a) cos” (1a) 
N, pa(L/a)*(1/h - 2°) cos 26 [2 = x/L] (10) 


--2 pa(L/e) sin 2g (1c) 


+ 


These expressions are derived under the assumption that the end 
stiffners at x = + L/2 cannot react longitudinally or torsional- 
ly. 

Calling: 


= longitudinal displacement, positive in the direction of 
ositive x 
2 See, for instance, ref. 2. p. 110. 
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tangential displacement, positive in the direction of 
postive ¢ 


radial displacement, positive inward 

shell thickness 

Young's modulus 
and letting Poisson's ratio VY = 0, as practically permissible 
for reinforced concrete shells, the displacements u, v, w are 


found to satisfy the following set of differential equations: 


(2a) 


(2b) 
(20) 


These equations may be successively integrated and lead to the 


following values of the displacements: 


2 
= - 22 (2)? cos - 27/3] (5a) 


2 
- 3 (Ey! sin (2*-1/16)- (27-1/4) 
(3b) 
+ 6 cos apt t(2"-1/16)- (2°-1/4) } 


2 


1 2 
(2)" gf (3e) 
under the assumption that: 
u(O0) =O; v(+ L/2) = 0 (4) 
It is seen from (1) that: a) Ng and Neg are zero along 
the boundaries ¢ = + n/2, i.e., that the horizontal edges are 
free of stress; b) N, is zero at the ends x = + L/2; 


3 See, for example, ref. 2, p. 117, eqs. (83). 
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c) the shear Ng at the ends x = + L/2 has the value: 


+ 1/2" + pa (L/a) sin 2¢ (5) 


Hence, the load p is carried by the membrane stresses (1), as 
for as equilibrium is concerned, provided the end stiffners 
develop the shears (5) (Fig. 2). 

On the other hand, membrane stresses will only develop in 
the shell if the boundary displacements (3) are consistent with 
the boundary conditions. Hence, membrane stresses will develop 
if and only if the ends of the shell: a) are free to move 
longitudinally; b) are prevented from moving tangentially; 
c) are allowed radial displacements and rotations; and if 


d) the horizontal edges of the shell are entirely free to move. 


36 
The spherical shell of Fig. 3 carries a uniform vertical 
load. We let: 
= meridional force per unit length, positive when tensile 


parallel force (hoop force) per unit length, positive when 
tensile 


uniform vertical load per unit of horizontal area 
rajius of srell 
colatitude 


longitude (measured from fixed axis coincident with 
cylinder axis) 


The membrane stresses in the spherical shell are given by: 


Ng ze 5 pa (6a) 


cos°¢) 


2 


while the shear is zero because of symmetry. 


i See ref. (1), eqs. (209), (210) with r, = 
2 
R= pa’ sin ¢, Z= pcos ¢. 
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Calling: 


meridiona] displacement, positive in the direction of 
positive g 


radial displacement, positive inward 
thickness of shell 
= Young's modulus 


letting Poisson's ratio Y= 0, the membrane displacements 


c 


given, by: 


(7a) 


2 
en ( 2 cos* ) (7b) 


while the parallel displacement u is zero because of symmetry. 
The disp?acements (7) are obtained under the assumption that 
v = at the hcrizontal edge = 

It is seen from (6) that a quarter-sphere (a sphere seg- 
mort bound by the diamctral planes $ = n/2, @ = n/2) will be 
in equilibrium under membrene stresses, provided; a) the 


reaction av the horizcntal edge be a vertical force: 


1 


b) the reaction at the vertical edge be a horizontal force; 


p2 (1-2 cos” g) (8b) 


iT 


5 = n/2 
(Fig. 4). 
Moreover, the quarter-sphere must be allowed edge displacements 
and rotations compatible with (7). 
lk, INTERACTION BETWEEN HALF-CYLINDER AND QUARTER-SPHERE 

The horizontal force (8b) on the vertical boundary of the 


quarter-sphere creates an equa] and opposite reaction in the 


cylinder; hence a longitudinal force; 


See ref, Pp. 370. 


E 
2 
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pa( 1-2 cos* g) (9) 
which is independent of x, must be added to the force (1b) in 
the cylinder due to the load (Fig. 5). 

Similarly, the shear (5) on the vertical boundary of the 
cylinder develops in the quarter-sphere a state of stress, which 
may be obtained from a solution for membrane stresses in 
spherical shells due to boundary loads. Flugge proved that? 


under forces: 


(10a) 


(10d) 


applied to its boundary 6 = n/2, a hemisphere develops 


membrane stresses: 


C n 
Ny =-Nj=- x (tan 6/2) cos n 


For n = 2, these equations reduce to: 


Ng = - Ny = - 299.2 (12a) 


cos. 6/2 


(12b) 


Hence, setting the boundary shear in the sphere (Nog of (12b) 
at 6 = x/2) equal to minus the boundary shear in the cylinder, 
(Ng at x = L/2 of (5) J, the constant C, takes the value 


6 See ref. 3, p. 43, or ref. 1, p. 379. Here the angles 


and @ are interchanged with the angles in the references. 
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pa (L/a), and the boundary shear (5) in the cylinder 


to develop in the quarter-sphere the forces: 
L, cos 2 
= pa (=) (13a) 
cos. 6/2 
(13>) 


(Fig. 6). 
These forces require for equilibrium an additional normal 


force on the vertical boundary of the quarter-sphere: 


** pa (2) cos 2g 


This additional normal force is equilibrated by an equal and 
opposite reaction on the cylinder boundary; hence a long- 


itudinal force; 


2) = - pa (2) cos 2 g (14) 


which is independent of x, must also be added to the force 
(1b) due to the load (Fig. 7). 


BENDING STRESSES AT BOUNDARY BETWEEN CYLINDER AND SPHERE 


At the boundary between cylinder and sphere; 
= 1/2 Nog Je = x/o 


he transverse and meridional forces Ng are unequal, 


L/2 (L/a) cos” [by (la)] 


= - (1/2) pa + (3/4) pa (L/n) cos 2 g [by (6a) 
and (13a)). 
Hence a difference in circumferential] strain 46; will be 
present at this boundary, of magnitude; 
Ae,~m 6 *x=1/2~ "6 


= (14 (5/2 -5 cos 
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The maximum end minimum values of Aé¢ occur at g = 5 and 
g = 0, respectively, and are given by: 
1 pa 3 ¢L = 
* 5 {1+ 5 (15a) 
-~-ipa 
bs -11 (150) 


The corresponding differential radial displacement *. between 
cylinder and sphere is equal to; 

woe ; (16) 
hence, bending moments and transverse shears will develop in 
the cylinder and the sphere capable of producing an equal and 
opposite Wo and of thus equalizing the values of Ng in the 


cylinder and in the sphere, 


An approximate evaluation of the bending stresses due to 
the differential strain dey is easily obtained by noticing that 
these stresses are locaj and that it is therefore permissible; 
a) to determine the cylinder deflections as if the cylinder 
were strained symmetrically about its axis; b) to determine 
the sphere deflections ac deflections of the cylinder tangent 
to the sphere at the boundary; and hence c) to evaluate 
the stresses in cylinder and sphere considering each long- 
{tudinal strip of cylinder and sphere as a beam on elastic 
foundation of spring constant Eh/a. 7 

If moreover we assume that the displaccment Wo will be 
equally compensated in the spnere and the cylinder, (i.e., 
that the value of Ny be the average of "y = 1/2 and 

Ng de no bending restraint will be developed at thre 
boundary, and the dispiacement w in the cylinder and the 


sphere must satisfy the conditions 
9 


aw , 
wo) at (17) 
7 See ref, 1, Pp. 389. 
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where x = O indicates now the boundary between cylinder and 
sphere (Fig. 8), 


The expression for such displacement is: 


w(x) = ef * cos 6 xX (18) 
and the ccrresponding moment is given by: 
2 
= -Du, A ef sing x (19) 
where: me 3 
2 Eh 


The maximum value of e f » sin p x occurs at 
Ax 0.8 and equals 0.3223. Hence by means of (15), (16) 


and (20) the maximum and minimum values of M, equal: 


M = - 0.0233 [1+ (2) Jpah (218) 


x “max 


Mlw 


M, 


x = 7-0-0233 [1 - (E) Jpah (21b) 


It is seen that these moments will be relatively small for the 


usual values of the ratio L/a. 


6. BENDING STRESSES AT HORIZONTAL EDGE OF QUARTER -SPHERE 


An approximate evaluation of the bending stresses at the 
horizontal edge of the quarter-sphere may be obtained by means 
of the same simplifying assumptions made in the preceeding 
section. 

A meridional strin of the spherical snell is approximated 
by a strip of the cylinder tangent to the shell at the boundary 
and the cylinder is assumed to be symmetrically strained. If 
the shell is built-in at the boundary a distribution of bending 


moments and transverse shears will be developed there capable 


B See ref. 
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ai 7 shells”, assuming that the law of variation of the displace- 
is | ments with @ be the same as the law of variation of the stresses, 
i.e., proportional to 1/ cos 6/2. 
Thus the meridional] displacement; 
o assuming v = O at g = 2/2, becomes by means of (13a): 
v= =) [sin 2 $+ sin 1n (tan ¢/2)] (24) 
salt as Similarly, the radial displacement is given by: 
aN, 
w= v cot 
(25) 
2 
= [6 cos" g-1+2 cos 1In(tan $/2)] 
cos 6/2 
The value of w at @ = n/2 is: 
a 


of producing radial displacements and rotations equal and 


opposite to those due to the loads. 
The radial displacement at the horizontal boundary due to 


the load p is, by (7b): 


2 


The boundary rotation is: 
1 


dw 
** 


(23) 


The displacement and rotation due to the stresses induced 


by the cylinder shears may be evaluated approximately by means 


of the displacement equations valid for symmetrically loaded 


9 See ref. 1, p. 370 
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and the value of the rotation is: 


1 dw 


** (27) 


Since d. is maximum at © = n/2, the maximum radial 
displacement at the boundary of the quarter-sphere is, by (22) 


and (26); 
1 


) (28) 


( 


while the rotation V is 
The displacement in a semi-infinite strip of cylinder 


zero. 


built-in at the end x = 0, i.e., with: 


w(o) = (29) 


is given by: 
- px 
w(x)=-de (cos 6 x + sin 4x) (30) 
and the corresponding moment is given by: 
2 _- x 
= 2D, Sef (cos @x - sin x) (31) 
where ¢ and D are defined py (20). 


The maximum value of * (cos gx - sing x) oceurs at 


x = 0 and equals 1; hence the maximum value of the bending 


moment, by (28) and (31), equals: 


re) 4 3 L 2 
M = - 0.144 [1+ 5 (=) pah (32) 


It is seen that this moment is approximately six times the 
moment developed at the boundary between cylinder and quarter- 


sphere, 


7. BENDING STRESSES AT HORIZONTAL EDGE OF BARREL 


The cylindrical barrel carries the load p by means of 


the membrane stresses (1), (9) and (14). Stresses (1) are due 


to the load p, while the last two stresses are due to the 
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interaction with the quarter-sphere and are independent of x. 
The displacements due to (1) are given by (3); the displace- 
ments due to (9) and (14) are obtained by integration of 
equations (2) with 

(2), (2) = pa [(1-cos" g)- (L/a) cos 2g) (33a) 


Ng = 0 (330) 


have the values: 


2 
5 BE (E) 2 [1 - 2 cos g - 2 (E) cos 2 g] (3ha) 


4 
2 


pa” Lye (Ly) 2° - 
(12 + 5 - 1/4] cos 2 g (34e) 
if conditions (4) are satisfied at the boundary. 

It is seen from (3) that the transverse and radial displace- 
ments, and the rotation due to the load p have at $= x/2 


and x = O the values: 


(35a) 


It is seen from (34) that transverse and radial displace- 


ments and the rotation due to y ) and x (2) have at 


¢ = x/2 and x = O the values: 
(362) 
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and 
x = 0 
- pe (by)? 2 (Lye 55 
4 4H, Eh (=) (5 + (=) (55b) 
x = O 


(36c) 


f the cylinder were supported by perfectly rigid end 
stiffners, the displacements at the middle cf the longitudinal] 
edges would be given by the sum of (35) and (36). The 
flexibility of the quarter-spheres introduces an additional 


vertical displacement MV_, which may be evaluated approxi- 


mately as follows: 


The crown of the cylinder ends moves down, due to the 
L 
(3c) at x = £ g = 0 by: 


v= (37a) 


The top of the sphere at the intersection with the cylinder 
moves down, due to the load and the cylinder shear reactions 


[eas. (To) and (25) J, 


2 
5 


The additional lowering in the cylinder is, therefore, 


(37c) 


The bending deflection at the intersection of cylinder and 


sphere introduces at ¢ = © 


{eqs. (15b) and (16) } 


1 


a downward displacement 


Thus the total downward displacement at the crown of the 


cylinder ends is: 
pan 


The vertical displacement at the middle of the horizontal 
edges differs from Wy, because, due to the difference in the 
values of Ng at the end and at the middle of the cylinder, the 
10 To obtain this value the limit of v cot ¢ mst be taken 

in (25) as $ approaches zero. This limit equals 


_ 
Fh ‘a 


| 
if 
| 
1 9 
cos 
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cylinder are from the crown to the edge shortens more at the 
midile point x = 0 than at the end x = L/2. The Ng at the end 
of the cylinder is now the average of Ng l, = 1/2 and 


2 


1 
Ng *2%6 x/2 


= pa (1 + ( 3/2 - cos’ g ) 


The corresponding strain is d = Ny/en and the transverse 


displacement at the edge becomes; 
n/2 pa” L 
af (1+2) (37f) 


a 


= L/ 
$= x/ 
The strain €y at x = O is equal to Ng/Eh, where Ng is 
given by (la), and the corresponding v is given by: 


(37g) 
n/2 
The difference between (37f) and (37g): 
9 
7, 
= 1/2 8 th (1 (=) ] (37h) 
= n/2 


must be added to w, to obtain 


due to the flexibility of tae supports: 


2 
+ ve = +f) +3 (BQ - } (38) 


If the horizontal edge of the barrel is built-in, a 
distrioution of statically indeterminate transverse reactions 
T, radial reactions R and moments M (all per unit length) is 
developed capable of inducing in the shell displacements and 
rotations equal and opposite to the sum of the boundary dis- 
placements and rotations (35), (36) and (38). T, R and M are 


assumed to vary according to a sine law in the x-direction. 
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The evaluation of displacements and rotations due to 
boundary forces in a cylindrical shell was until very recently 
@ non-elementary and very tedious task requiring the integrat- 
ican of an eigth order differential equation; bu*t thanks to the 
jabors of Mr. A. Parme it is now possible to evaluate such 
quantities in a routine manner by using the tables of the ASCE 
Manual of Engineering Practice No, 31, "Design of Cylindrical 
Concrete Shell Roofs" (ref. 2). 

Thus the Parme manual contains tables for the direct 
evaluation of the displacements and rotations at the middle of 
the cylinder edge due to a transverse unit force 
T = lxsin x x/L, which we indicate by Vipr Wp and @,, and the 
same displacements and rotation due tc a unit radial force and 
a unit moment, which we indicate by the symbols Vp» “pe st and 


Ye and respectively. 


Equating the displacements due to the loads to the dis- 
Placements due tc the statically indeterminate reactions of 
amplitude T, R and M (as is done in any statically indeterminate 


structure) we obtain the set of three simultaneous equations; 2+ 


Vn T + YR R + YM M AV. (39a) 


+, AH, Jo) 


6, T+6,R+6,™ (39c) 


which allows the evaluation of the reaction amplitudes T, R 


and M, An additional table of the Parme manual permits then to 


determine directly the maximum value of the bending moment due 


11 The positive sign of v and w in this paper is opposite 
to the sign in the Parme manual. The positive sign of 
® is the same, 


| 
293-15 = 
| 


8. EXTENSIONS OF THE PRESENT THEORY 


to these reactions, which in most cases equals i 


In this manner, although explicit formulas for them can- 
not be given, the bending stresses at the horizontal edge of 
the cylinder may be determined in an easy fashion and the 
analysis of the composite shell may be completed without great- 
er difficulties than those entailed by the analysis of any other 


statically indeterminate structure, 


The methods of this presentation may be extended in a 


variety of ways. 


a) The analysis of the same type of composite shell under dead 
load follows step by step the present analysis, and formulas 
analogous to those derived here may be easily written, witn one 
exception, In the study of the stresses in the quarter-sphere 
duc to the cylinder shears it is found that the shears vary as 
the sin ¢ and that in (10) and (11) the integer n equals 1. 
Hence at @ = O the stresses defined by (11) become infinite, 
due to the fact that the shear distribution is not equilibrated 
over the complete circle. This difficulty is avoided by 
considering an equivalent equilibrated shear distribution, 
obtained by expansion of sin gp into @ sine series containing 


only even sines; 


: 
sing=+ (-1)"* sin 2m (ho) 
m=] 4m -1 


and by using the first few terms of this series in tne deter- 


mination of the corresponding stresses and displacements. 


12 It is possible to refine this treatment of boundary stress- 


c 


es by taking into account the first two terms of the sine 
series expansions of the reactions T, R and M. Additional 
tables in Parme's manun] allow this refinement, but quite 
often the first term of the series will suffice to indicate 
the order of magnitude of the maximum bending moment in the 
vicinity of the horizontal edge. 
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b) The same analysis with minor changes may be extended to 


segments of cylinders and spheres limited by an angle ¢ 


max 
less than x/2. Hence composite shells with less than x/2 


may also be considered. 


c) Loads uniform across the shell or varying transversely as 
the dead load, but varying along the shell length may be taken 


into account by means of Fourier series expansions. 


d) Stresses due to uniform temperature changes may be deter- 
mined by substituting for the displacements due to the load, 
the change in lateral dimensions due to temperature changes. 
Thus, for example, the bending moment in the quarter-sphere 
duc to a change of T °F is obtained by substituting in (31) 


for the value: 


where o¢ is the coefficient of thermal expansion. 
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